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Abstract 

The connection between symmetries and linearizations of discrete-time dynamical systems is being inverstigated. 
It is shown, that existence of semigroup structures related to the vector field and having linear representations 
enables reduction of linearization problem to a system of first order partial differential equations. By means of 
inverse of the Poincare map one can relate symmetries in such linearizable systems to continuous and discrete ones 
of the corresponding differential equations. 



1 Introduction 

The principal objective behind the concept of (exact) lin- 
earization is to express an apparently involved evolu- 
tion law in simple terms, facilitating further analysis. 
The purpose is also to bring a system to some 'canoni- 
cal' form enabling comparison among studied dynami- 
cal systems in terms of e.g. minimal embedding dimen- 
sion, which can be regarded as a complexity measure. 
To that end one seeks a projection map p : $ — * M from 
a suitable linear space S onto the phase space M, which 
defines a homomorphism of appropriate (semi)group ac- 
tions. The desired feature is of course finiteness of the 
embedding space S, otherwise the linearization remains 
formal. On the other hand, it is worthwhile to recall 
the applications of linearization in infinite-dimensional 
spaces. An example are the methods of Hilbert spaces in 
spectral theory 1 1 [. 

Since the works of Lie, the notion of a symmetry as 
an automorphism in the space of solutions plays cen- 
tral role in the analysis of differential equations |<3|7||H1. 
In particular, it is well-known that the knowledge of an 
infinitesimal Lie-point and generalized symmetriy, al- 
lows to reduce the system of ordinary dyfferential equa- 
tions (ODEs). An analogous result for ordinary differ- 
ence equations (OAEs) was obtained by Maeda |2[, and 
expanded by Quispel etal. 1 3 4 5 1 . Namely, they showed 
that determination of a Lie-point or generalized symme- 
try, enables reduction of the dimension of an OAE, and 
in particular cases its linearization |2 4J. The approach 
taken up by Maeda and Quispel et at relies on finding 
an infinitesimal time-translation symmetry, and trans- 
formation into normal coordinates. It should be noted, 
that determination of such particular symmetry may be 
in general highly nontrivial task. 

In this work, we introduce a method for the study of 
difference equations, based upon observation that finite 
symmetries of the evolution generator form a structured 
set of relations, identified as a semigroup of similarities. 
If at least one of its closed subsets is known, then that 
structure can be used to form a linear semigroup of iso- 
morphic structure, and obtain differential system of con- 
straints for the desired projection map p. This method is 
illustrated by two examples: the 'classical' Riccati equa- 
tion related to the symmetry group SL(2), and a thor- 



ough analysis of a logistic equation, with emphasis on 
the role played by similarities in classification of ob- 
tained continuous-time solutions (the inverse Poincare 
maps). 

2 Linearization and symmetries 

Let us begin with a brief account of the problem of lin- 
earization of recurrences, and its relation to the symme- 
tries of evolution generators. 

Let / : ./# — > ^# be a function defined on a mani- 
fold ./#. Consider the autonomous system of difference 
equations 



x n +i=f(x n ), neN o = NU{0}. 



(1) 



Iteration of the function / generates an associated semi- 
group with multiplication defined by composition of 
functions, and unity (monoid) 



G := {/" | / u = id, / 



n+l 



f n o f, n e No}. 



The linearization of Eq. Q is a pair (L,p) comprising a 
linear map L : § — > $ acting on a vector space S ', and 
projection p : S —* y£ , such that the following diagram 
commutes: 



S 



S 



r 



for all n G No. Equivalently, one can write 

f n op = poL n . (2) 
Iterations of the linear map L form another semigroup 

G Un := {L n | n E N }, 

which, acording to 10 is homomorphic to G. Therefore, 
the linearization ivolves (i) fibration § — > jf( , (ii) linear 
function L, and (Hi) semigroup homomorphism. 

Further, the kernel of p is a free monoid acting along 
the fibre: 



kcrp = \k: $ — > (o \ p o k = p} 



(3) 
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By definition, the elements of this kernel satisfy 

V K 3 K ' : k o L = L o k (4) 

where k, k' G kerp. Notice, that the sections {a} £p _1 
and consequently the kernel elements too, are typically 
multivalued. Whenever they appear in some expression, 
it is understood that the equality holds true only for par- 
ticular subset of these maps. For instance, the inverse to 
Eq.0 

a o j = no L o a, 

means that for every a G p^ 1 there exists k G kerp (or 
vice versa) satisfying this relation. 

The iteration of function / is often a result of tempo- 
ral discretization, and interpreted as the Poincare map of 
some underlying, continuous-time system. The recipro- 
cal operation is what we call the inverse Poincare map 1 
(IPM), namely an extension (/ n ,^#) (/*,^) ; such 
that Jt D J£, /* : ./# — » for t G K+, and 

(0 f 1 = f, 
(ii) ]im t \o f = id, 
(m) VMel+: J t o/ s = / t+s , 

In other words, the inverse Poincare map is a one- 
parameter continuous monoid {/*} congruent to {/"} at 
t G No- The linearization (L, p) of Eq. provides a suit- 
able framework for obtaining IPMs. It must, however, 
be borne in mind that unlike L n , the power 1} , t G K+ 
is multivalued what makes characterization of the semi- 
group Gun by its generator L imprecise. One should in- 
stead use a particular branch of its logarithm E, that is 



L* = exp(tE), 



which is well-defined. If the kernel of p is linear, then {3} 
implies n o E = E o n'. 

Further, the notion of inverse Poincare map allows 
to discriminate between apparently different lineariza- 
tions, namely two linearizations (pi, Ei), (p 2 , E2) are re- 
garded as equivalent if the corresponding IPMs 

ft = Pi ° e' Bl o ai, 
fi = P2<>e tE *oa 2 , 
are homomorphic, i.e. 

3 4> : JZi =$ .M 2 : 4> o f{ ^ /* o 0, 

wwere the symbol indicates multivalence, and '=' 
is used instead of '=', to emphasize that the branches of 
map </> on both sides may be different. 

According to the above definition, equivalence im- 
plies that, besides <f>, there must exist a map U : S\ — > S 2l 
such that 



p 2 O U = (popi, 

k 2 o e tE2 o[/=[/o e tEl 



(5a) 
(5b) 



For linear kernels the condition 15bt reduces to 'similar- 
ity' between the logarithms k 2 o E 2 o U = U o E^ a K\. 



Whenever the two flows are equal, f{ 
of eqs. is of the form 

U cq — <y 2 Pi , 



/J, the solution 



(6) 



and <fi = id. Indeed, this reduces eq. t5a> to an identity, 
while i5bl can be verified by direct calculus: 



K 2 oe 



~tE 2 



oU ei 



— <J 2 ° p 2 ° e 

= ^2 O f\ O Pi 
= (T 2 o ft o px 

= G% ° Pi ° Ki O C 



o <7 2 Pi 



'-Ei 



t/e, 



If /* 7^ /!/ an d one can find an alternative solution of 
<5B . that is 17 ^ t/ cq , then 



p 2 oU O (Tl , 



(7) 



is one of the desired homomorphism maps, establishing 
equivalence between these two IPMs. 

At the same time <f>, becomes a symmetry of the func- 
tion / itself, or of its extension defined on ^#2, that is 



' / = / ' 



(8) 



Conversely, any such discrete symmetry provides sim- 
ilarity transformation between two IPMs. If the image 
of <f) is j$ , then it is a proper symmetry of the original 
function /. Equation directly suggests another way of 
finding symmetries by means of the linearization (p, L). 
Any linear mapping M : § — > S whose commutant with 
L is an element of the kernel of projection map 



La M oT 1 0M- 1 e kerp, 
gives rise to a symmetry <p defined by 

4> := p o M o a. 



(9) 



(10) 



The case, M — L l for fixed t G R+ leads to the 'triv- 
ial' symmetry coinciding with the time-translation one 

<P = f- 

The condition JSJ, can often be rewritten in alterna- 
tive form as a solution of certain implicit system 

S(x,<f>(x))=0, 

such that the following implication holds 



S(x,x') = 



S(f(x),f(x')) = 0. 



(11) 



1 We use this term, since it ascribes certain physical meaning to the definition, however from mathematical viewpoint, a more apt name would 
be 'monoidal homotopy' - homotopy /' between /" and f n+1 which additionally has the structure of a monoid. 
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Accordingly, the solutions of fj) are invariant subsets of 
this relation S(x n ,x n+ i) = 0, i.e. the symmetry maps 
solutions into solutions. 

The symmetries 4>{x), similarly to linearizing trans- 
formations p, can be expanded locally into power series 
in x in the neighbourhood of fixed points (FPs) of the 
generating map /. Moreover, the same reasons which 
cause nonexistence of (local) linear normal forms, also 
prevent unambiguous symmetry expansion whenever / 
has resonant spectra at at FP x*. If we denote by x % , 
i = 1,..., dirru/# the local coordinates of a point x G . // , 
and the spectra of the first derivatives p.j = dpjdx^ by 

sp f\j(x*) = {A(i),..., A (m) }, to = dim^, 

then this resonance condition takes the form 

dim ^ 

BfceN™ : [] (A (i) ) fc * espf, 
1=1 

If the condition JS} is to be satified by a single branch of 
<\> near FP x*, then one finds <fi(x*) = x*. The expansion 
coefficients in the series 

<f>(x* + Sx) 1 = xt + (j)\ 3 8x J + \(j)\ 3k 5x j 5x k + ..., (12) 

where summation over repeated indices is understood, 
are recursively defined functions of f\j (x* ) and the free 
parameter tensor (x*). In particular, chosing 



4>*A X *) = (f n YA x *)> 



n G Z, 



leads to the trivial symmetries <fi = /", n e No. Never- 
theless this does not imply that (f> coincides with the time 
translation symmetry /' for all t G M + , since in general 
x* is not a fixed point of /*. Other special choices in- 
clude </>*.,- (x*) =0 => <j) = x* = const, 4> l j(x*) — S l j 
4> = id, and </> l j(x*) = —S i j yields self -inverse symme- 
tries, i.e. <\>~ x = (f>. In the second example provided in 
this article, we use local expansion <12> , together with the 
property dlli . to investigate numerically the symmetries 
of logistic map. In the same example, the many- valued 
symmetries do not satisfy 4>{x*) — x* for all fixed-points 

Composition of symmetries form semigroups with 
identity (groups, if all elements are invertible), which 
can be in general uncountable. The whole set of symme- 
tries of a given map / is always non-empty, containing 
at least the semigroup G itself. In particular, if a subset 
5^ := {4> s } admits continuous parametrization s G R, 
satisfying lim s ^o 4> s = id, 4> s ° 4> l — <fi s+t , then one can 
construct a piece-wise autonomous (i.e. autonomous on 
open intervals of the evolution parameter) vector field 



s=0 



<P(x), 



$ := ¥d t + <t l d; 



where barred indices refer to complex-conjugated vari- 
ables if ^ is a complex manifold - real case is obtained 
by restriction to vanishing imaginary parts. Differentia- 
tion of the symmetry condition 



leads to the tangent map (push-forward) @\f = f*<P, 
which in local coordinates has the form 



(13) 



provided / is holomorphic. As pointed out by Maeda 
|2J and Quispel et ah (3] [3, existence of the vector field 
<fi satisfying ( 1131 . can be used to reduce dimensionality 
of Eq. £Q in the same way as it is done in the theory of 
differential equations. More precisely, if there exists an 
IPM /' such that the symmetry condition holds true for 
all t G R+, then the corresponding vector field X defined 
by 



dt 



X{x), X :=X l d t + X l d- l 



satisfies the Lie point-symmetry condition 



[X,$] = 0. 



3 One-parameter family of similarities 

Instead of constructing symmetries and similarities 
upon linearizations, we now consider a situation when 
one knows at least some of them in advance. The ques- 
tion arises as to whether this information can be ex- 
ploited to determine the linearizations, and in this para- 
graph we investigate this possibility. 

The generating function / together with its sym- 
metries form a semigroup with respect to composition, 
equipped with unity. That the function in concern / may 
be part of a structured set of maps, is a key element, 
which significantly simplifies the search for a homomor- 
phism with linear representation of the same structure, 
that is - a linearization. 

Similarity between two maps generalizes the notion 
of a symmetry. Analogously, a set of similarities form a 
semigroup with respect to composition 



H := {h(a) : Jt =4 Jt \ 

V a!1 a 3 7 : h(a) o h(j3) 



h(j)oh(a)}. 



In particular, if all elements of H are inverible then this 
set becomes an inner automorphism group with compo- 
sition h(a) : h{/3) h-> h(a) o h(J3) ° h~ 1 (a) = h(j). 

Suppose, that H is continuously parametrized, and 
contains the function /. The structure of this set can 
be used to determine (formally) isomorphic semigroup 
Hun = {L(a}} of linear maps acting on some vector 
space £ . Assuming there exists a homomorphism be- 
tween H and Hi; n the problem of linearization reduces to 
finding a map p ; $ — ► j% such that 



'/ = /< 



V Q : p o L(a) = h(a) o p. 



(14) 
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Let the explicit, coordinate form of L{a) be 



L a (y;a) = L\(a)y b , 



a, b = 1, 



,m, 



L a b {a) e C. 



Differentiating with respect to y and a gives two inde- 
pendent relations 

p^Ly^^KMvWAv) 

p\ a (Ly)d a L\y b - d a h\p{y)) 

where y = {y a } is the coordinate vector of the point y, 
d a = d/da. Comparing the above equations we get 



h\ ]P >, b d a E b c y c = d a h\ 



(15) 



where all omitted argumets of p are y, and as before 
L( a ) = e E( - a \ that is d a E b c = [L~ 1 ] b d d a L d c . There- 
fore, the existence of continuous 1 -parameter semigroup 
of similarities leads to the first order system of partial 
differential equations (PDEs) for the mapping p. Its solu- 
tion, together with the linear map L(a) such that h(a) = 
f, provides the desired linearization of the system 

We stress, that the crucial point in the above construc- 
tion is the formal isomorphism between structured sets 
of nonlinear and linear mappings, and not their particu- 
lar composition rule - cf. example B . 1 I below. 

3.1 Example: homographic map 

For an easy illustration of the procedure described 
above, we first analyse the case of Riccati equation. Let 
x G U = C, and the dynamical system fl} be the equation 
generated by a homography / : U — > U 



f(Xn) ■-- 



ax n + (3 
7X„ + 5 ' 



Recall, that the set of all homographies is a Kleinian 
group, and under the above constraint imposed on 
the parameters, which exclude uninteresting constant 
case, it becomes isomorphic with SL(2,C). A local 
parametrization can be given by 



L 



j(8 + jx) 

ctS-l 



6^0 



(17a) 
(17b) 



Clearly / G H. 

Out of the three continuous parameters (a, 7, S) it is 
sufficient to take the one specific - choose a for concrete- 
ness. The differential constraint fTpl becomes 



dh dp 



d 



dh 

oa oa 
Inserting llT7l . one gets single PDE of the form 

dp S dp 7 1 



dy 1 (6 + 7p) 2 dy 2 (S + jp) 2 jy 1 + 5y 2 



The solution of this equation has two integration con- 
stants and D: 



p(y) 



SCD 



(18) 



Taking into account ITU , IT71 and fT8l , one finds that the 
only possibility left to satisfy C 5^ is to set D = 0. 
Therefore 



p(y) 



which is the desired linearizing map. 

Having foimd the linearization, it is straightforward 
to provide the inverse Poincare maps /* (x) of the Riccati 
equation 116\ . These maps are also homographies with 
t-dependent coefficients, and in general, for arbitrary ini- 
tial position x G U they go to infinity in finite time. This 
implies, that the motion is better described in complex 
projective space jjt = CU{oo} = CP 1 , instead of a sin- 
gle chart U. Taking the second local chart as U' = C, 
with transition functions g : U — > U' x' = g(x) := 1/x, 
one has /' = g o f o g^ 1 , that is 



-1 = no = 



l(5x' n +7) 
(aS — l)x' n 



The projection map in in these coordinates p 1 = g o p, is 
simply p' (y) = y 2 /y 1 . Both p and p' define local inhomo- 
geneous coordinates of the projective space J£ = CP 1 . 



(16) 3.2 Example: the logistic equation 



Consider a more involved example of the logistic equa- 
tion 



^n+l f\Xn) • 4x n (l Xyi). 



(19) 



It is well-known that it is linearizable - we now recover 
this linearization following the approach described ear- 
lier. We first observe, that / is a member of the following 
set 2 



(20) 



2 The branches of infinitely multivalued map arccos : [—1,1] =f M, 
- l)/2], where Arccos : [—1,1] — » [0, tt] is the principal branch. 



H = {h(x; a) = sin 2 [77a + arccos(l — 2x)] \ a G M} 

for aeZ. The compositions 

h{a)oh{(3) ^/i(7)o/i(a), 

satisfy either 

7 = 2/3 — a or 7 = 3a — 2j3. 
can be expressed explicitly as arccosfc(ic) = (— l) fc Arccos(x) + Tr[k 
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The following linear operators L± = L±(a) 



L+ = 



and L- = 



define two representations Hii„ of 1201 . that is 

L+(a)L+{0) = L+(2/3 - a)L+(a), 
L_(a)L_(/3) = L_(3a - 2/3)Z,_(a). 

The equations l !15> . in each case reduce to a single PDE 

dp . 2ir 



dy 1 



where the '±' signs correspond to i+ and L_ respec- 
tively. Their solutions however, are formally identical 



p(y) 



C± = C±(y 2 



Substituting p(y) into Eq. iTH p{L±y) — h(p(x)), we find 
the functions C± to be constants, and for each label '±', 
two further possibilities arise, namely we have 



p+Av) 



\ cos (tt(^ 



P+,a(y) = 3 - !cos(tt( 



-«)). 
-3a)), 

F-,i(2/) = |-| cos ( 7r (^ 

-a)). 



p-Av) = \ 



5 COS (7T(^ 



The corresponding inverses are 

°'+,i( a; ) = [y 2 (^arccos(l - 2x) + a),y 2 ], 

0-+,2(x) = [y 2 (iarccos(l - 2x) -3a), y 2 ], 

o"_, 1 (ar)= [2/ 2 (iarccos(l - 2x) + ±a),y 2 ], 

o_ 2 (o;) = [2/ 2 (iarccos(l - 2x) -a),y 2 ]. 

The nonlinear IPMs can be computed from 

f ± = p o L± o a 
where the powers of L± are 3 

a(2* - 1) 



(21) 



L* = 



2' 




1 



(-2)' f(l-(-2)<) 
1 



Now, the composition J21> yields 

f +1 (x) = sin 2 (2 t_1 arccos(l -2x) +na{2 t - 1)), 
f +2 (x) = sin 2 (2'" 1 arccos(l - 2x) + na(l - 2*)), 

f-A x ) = sin2 ((-2)* -1 arccos(l - 2x) 
+fa(l-(-2) t )), 



/' 2 (*) 



sin 



((-2)*~ 1 arccos(l - 2a;)]). 



Clearly, /l and /I 2 represent the same family of maps, 
and for the second pair, /i 2 is a subset of ft i- There- 
fore, there are actually two distinct sets of (multival- 
ued) continuous trajectories ft and ft, however only in 
the first one there exists a real-valued branch of 2*, and 
hence f\{x) G K, provided x G [0, 1]. All remaining tra- 
jectories are complex. 

Our purpose now, is to discuss the symmetries con- 
structed from obtained linearizations. Since £±,i = 
L±,2, therefore we use U = 1 as a solution to i5b\ . The 
first pair composed according to Eq. that is 



P+,i o 0+.2 = p+,2 o cr+,i 

• [ira + h arccos(l — 2x)] , 



= sin 



becomes an identity for a G Z, which is a trivial symme- 
try. But the second pair 

° = P-,2 c-.i 

= sin 2 [§7ra + ^ arccos(l — 2a;)] 
=: <p e (x) 



yield a nontrivial one at a 
duction this is 



(f> c ( x ) = l(z-2x± 2y/3x(l-x)j 



1,2 + 3k, k G Z. After re- 
(22a) 



namely, the solution a;' = O (a;) to an equation of ellipse 
inscribed into the square [0, 1] x [0,1]: 



S e (x,x') := (x + x' - l) 2 + \(x 



x'f -4 = 0. (22b) 



Compositions of 4> e with /± gives uncountable number 
of symmetries out of this one. Nevertheless, this set is 
not exchaustive; other symmetries can for instance be 
constructed according to Eq. I llOt . In order to illustrate, 
consider the following linear maps 



AL 



+ - 



i 


;i - l)a 


M_ := 


"i |(l-i)a" 





1 




1 



trivially satisfying the condition by [M±,L±] = 0. 
Restricting to the case of integer parameter a, the com- 
positions p± o M± o a± =: 0±.h take the form 



</>+,h = — sinh 2 \ira + 2 arccos(l — 2a;)] , 

c/>_ h = - sinh 2 [^(1 + i)a + \ arccos(l - 2a;)] . 



We point out that, neither of these map the interval [0, 1] 
into itself, as <p e does, therefore they are not proper sym- 
metries of the logistic equation in the sense of definition 
given earlier. In particular, both the above sets have one 
common element for a = 0. It is rather special one for 
being real-valued and self -inverse (f>7 1 = 



3 In the complex domain both 2 t and (—2)' are countably many-valued. A particular branch E± 

in (±2)' = exp[t ln2 + ti7r(2fe + (1 =F l)/2)]. 



In L± is selected by fixing the integer k 
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(sin 2 [i arccosh(l — 2x)l 
inh 2 [-1 arccos(l — 2x)] 

_l„2 n /n_ -, \1 



- sinh 2 ^2 ■ 

cosh 2 [i arccos(2a; — 1)] 



[cos 



arccos(_ 
i ^2 arccos(2a; 
arccosh(2a; - 



x e (■ 
*e[o,i], 

l£ [1, 00 



00,0], 



This relation 0h, together with <f) c , are illustrated in 

pig. mo- 



u = 4 





4>e 


0h 


f 


\ 





a) 



1 

a; 



/i = 3.83187 




b) 



/i = 3.5625 





0c \fl 




f 







1 

.7' 



/i = 2 





\0h 


f / 

J / 


\ \ 



d) 



Figure 1 Examples of the symmetry curves of general logistic 
equation f(x) = lix(1 — x), in several cases of its control parameter p. 
Due to / o s = f, s(x) = 1 - x, both </> c and </>h have point-reflection 
isometry through (a:, a:') = (i, |). The two cases n = 3.83187 (close 
to periodic window) and ju = 3.5625 (just before onset of chaos), were 
found numerically through procedure described in the text. The curves 
4>c on b), c) seem to be fractal, for not being closed and filling densely 
subsets of ..iff 2 = R 2 . 

We have also investigated symmetries of the logis- 
tic equation x n+ i = /J,x n (~l — x n ) numerically, in other 
cases of the parameter fj,. The applied procedure directly 
uses the symmetry definition ( II li . Namely, assuming 
</>(x») = x*, one can develop into power series iT2i in 
some neighbourhood of the fixed point described 
in Sect. |3 We confine ourselves to the self -inverse sym- 
metries _1 = 0, by fixing the free expansion param- 
eter to d x <j)(x*) = —1. Starting from x := x* + Sx, 
\8x\ <C 1, the expansion can then be used to approximate 
x' := <f>(xo) to an arbitrary precision. By iteration 

{x n ,x' n } i-> {x n+1 ,x' n+1 } = {f(x n ),f(x' n )}, 

a finite number of points belonging to the symmetry re- 
lation is obtained. Many repetitions of this algorithm, 
with various initial displacements Sx yielded sufficient 
data to draw reliable pictures of in number of cases 
fi e [2,4]- two particular are shown on Fig.nj),c) 

Both O and 0h exist for all checked values of the con- 
trol parameter /i, except \i = 2 where C degenerates to a 
single point at (x*, x*) = (\,\)- At this extreme the sec- 
ond curve 0h can be given closed analytical form, thanks 
to the well-known linearization 

2x{l - x) =p(2a(x)), 

p(y) = l(l-e 2 y), a{x) = \\n{l-2x). 

Again, taking M = — 1 in Eq. JTOi gives (cf. Fig.QJi)) 



(f>h(x) = p(-a(x)) 



2x -1' 



(23a) 



or, in the equivalent form (cf. Eqs. Il22at . <22bt ). the solu- 
tion x' = 4>h(x) to 



Sh(x, x') = xx' — i(x + x') = 0. 



(23b) 
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It should be noted that, in contrast to all other cases of 
H > 2, here 4> h is an ordinary (i.e. single-valued) func- 
tion. 

We also remark, that by construction each of the fixed 
points x* = 0,1 - iT 1 of the vector field / is also a fixed 
point of one of the two symmetries 0h and 4> c respec- 
tively: 

O^img0 o (O), 1 - E img0 o (l - fi^ 1 ), 

6 img0 h (O), 1 - £ img0 h (l - 

This is not generally the case for the time-translation 
symmetries /±- It is therefore plausible to conjecture 
that, apart from special degenerate cases, presence of 
non-trivial discrete symmetries associated with equilib- 
rium points of the vector field / : .# — > ./# is a generic 
feature in difference equations. In this sense, <j> e and 
</>h can be considered typical of logistic map f(x) — 
[ix(l — x). 



[5] Sahadevan R., Byrnes G.B., Quispel G.R.W., "Linearization 
of difference equations using factorisable Lie symmetries", in 
Symmetries and integrability of difference equations ed. D. Levi, 
L. Vinet, P. Winternitz, CRM Proc & Lecture Notes 8 AMS, 
1996, pp 337-343; 

[6] Olver P.J., "Applications of Lie groups to differential equa- 
tions", Springer- Verlag, New York, 1996. 

[7] Bluman G.W., Kumei S., "Symmetries and differential equa- 
tions", Springer- Verlag, New York, 2nd. ed. 1989. 

[8] Stephani H., "Differential equations. Their solution using 
symmetries", Cambridge Univ. Press, Cambridge 1989. 



4 Conclusions 

In this paper, we have investigated the relationships 
among symmetries, similarities, linearizations and in- 
verse Poincare maps of ordinary difference equations. 
The general conclusion is, that the existence of addi- 
tional, continuously parametrized family of similarities, 
admissible by the system, plays an important role for the 
analytic solution of OAE and /or its linearization. Ex- 
istence of such structure can be regarded as an analog 
of Lie algebra in the case of Lie-type ordinary differen- 
tial systems. This structure generalizes the notion of a 
symmetry for it does not necessarily have to be a group, 
but may instead form a semigroup without unity. The 
important point is, that if a linear representation of this 
structure exists, then the linearization problem of origi- 
nal OAE can be reduced to a solution of a finite partial 
differential system. 

Finally, as suggested by the results presented here, 
the number of non-trivial discrete symmetries appears 
to be directly related to the cardinality of the fixed subset 
of phase space Ji . Interestingly, these relations need not 
to be simultaneously symmetries of the inverse Poincare 
maps /*, but instead provide similarity transformations 
among them. 
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